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Ps| ' We discuss the cosmological consequences of QCD phase transition(s) on the early uni- 

^1 verse. We argue that our recent knowledge about the transport properties of quark-gluon 

^ ' plasma (QGP) should thraw additional lights on the actual time evolution of our universe. 

^0 ! Understanding the nature of QCD phase transition(s), which can be studied in lattice gauge 

' theory and verified in heavy ion experiments, provides an explanation for cosmological phe- 

fS| I nomenon stem from early universe. 

1 introduction 

Our knowledge about the early and present universe is based on two successful models; the 
r~| \ standard model for elementary particles and the one for cosmology. The main ideas about 

the early universe cosmology seem to be confirmed from various cosmological data, like the 
temperature anisotropy in cosmic microwave background radiation (CMB), the light element 
\ abundances, etc. Although we have a solid description for phase transitions in almost all epochs 

' in the early universe, their confirmations still indirect processes. The observation of cosmological 

00 ■ 'exotic' phenomenon and the fruitful results from heavy ion experiments are promising tools to 

\ confirm these theoretical predictions. 

On I In this work, we limit our the discussion to QCD epoch, i.e. to energy scale of /S.qcd ~ 

. onn \^o,\r + ~ in-5 _ in-6 

00 
O 



200 MeV, t ~ 10~ — 10~ second after the begin and when the universe has the Hubble radius 
of about dff 10 Km corresponding to scales of about 1 pc or about 3 light years today. 
^ I At high energy the coupling between quarks and gluons likely becomes weaker. This coupling 

diverges at Aqcd = 200 MeV. From this property, QCD is known as asymptotically free theory. 
As the energy decreases, i.e. the universe expands and cools down, the quarks and gluons 
I underwent to confined state, to hadrons (deconfinement-confinement phase transition) and the 

chiral symmetry almost simultaneously breaks down, i.e., orientation of right- and left-handed 
quarks. The latter likely leads to production of Goldstone bosons. The three pions are the 
lightest Goldstone bosons. Dynamical fiuctuations on these particles and/or disoriented chiral 
condensates (DCC) would manifest this transition. 

In this talk, we discuss the cosmological consequences of nature and order of the QCD phase 
transition(s) on early universe. 



2 Ideal and non- ideal Quark Gluon Plasma 

A theoretical framework for dynamics of phase transition(s) from quark-gluon plasma (QGP) 
to hadrons still fails. In heavy ion collisions, we can detect produced particles at final state, 
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i.e. after chemical and thermal freeze-out. Consequently, there is no smoking-gun signal for 
occurrence of equilibration processes while the energy density e decreases. Same statement is 
also valid for the degrees of freedom (d.o.f.). It is known that the thermal (statistical) models, 
like hadron resonance gas, work very well in final state [BE]- These models presume a charge- 
conserved hadronic phase, but they have no access to the phase transition itself. 

Although the energy density e available to heavy ion experiments at CERN-SPS [3j, BNL- 
RHIC 15 and recently to CERN-LHC [5j exceeds the critical value calculated in lattice QCD 
(ec ~ 2 GeV/fm~^ for physical quark masses and vanishing net baryon number density), there 
is no unambiguous QGP-signature. For example, e achieved at RHIC ~ 5.5 GeV/fm'^ for proper 
time To ~ 1 fm/c. If Tc ~ Tjg, where Tc is the critical temperature and Tj^ is the freeze out 
temperature, holds at small or vanishing net baryon number density, as we used to assume for 
ideal QGP, the phenomenological signals characterizing the phase transition(s) have to remain 
measurable at final state, at least the ones which are not sensitive to strong interaction or to 
the medium, such as photons, leptons or color screening as J/^ dissociation into two leptons. 
According to recent lattice QCD results with almost physical quark masses, Tc ~ 200 MeV [6]. 
This value exceeds Tf^ with at least ~ 30 MeV. During the relaxation time, within which 
the system (universe) loses ~ 30 MeV most of dynamical fluctuations might be moderated. 
According to [7|, the inelastic scattering rate should have steep T-dependence. 

The plasma, like QGP, is a state in which the charges are screened due to the existence of 
other mobile charges. This will modify the Coulomb's law. Lattice QCD results in that the J/^ 
bound state would survive up to 2Tc, where Tc ~ 270 MeV [8]. The RHIC results suggested to 
include viscous corrections to the hydrodynamic evolution [9j , so that the plasma state turns to 
be a fluid rather than an ideal gas of massless components (quarks and gluons) and negligible 
correlations (free interactions). 

Taking into consideration shear viscosity in QGP fluid leads to slower hydrodynamic evolu- 
tion relative to QGP with zero viscosity. The transverse expansion in viscous QGP fluid turns to 
be stronger than the one in an ideal QGP fluid. Also the particle production might considerably 
be enhanced in viscous QGP fluid. That QGP turns to be viscous fluid, apparently, modifles 
our ideas even about the time evolution of early universe. 

3 QCD Phase Transition(s) and Chemical Freeze out 

The QCD predicts that the asymptotically free quarks and gluons are weakly correlated. QCD 
has been extensively studied on lattice for the last thirty years; QCD Lagrangian has to be 
discretized. Then we put everything on a finite a space-time lattice. It has been shown that 
a rapid change in various thermodynamic quantities undoubtedly takes place at sufficiently 
high energies [10]. The degrees of freedom markedly decreases in a relative narrow region of 
temperatures. 

First order phase transition at Tc ~ 270 MeV is evidenced in lattice QCD without dynamical 
quarks (quenched) [111 I12j . Including the dynamical quarks makes the order of phase transition 
and the value of Tc depending on the mass and number of quark flavors. 

For two massless light quarks (m^ = = 0) and inflnity heavy strange quark (mg oo) 
at vanishing net baryon number (chemical potential ^iq), which is appropriate to the early 
universe, where /iq/T ~ 10~^, the phase transition is of second order and Tc ~ 175 MeV. If the 
light quarks get small masses, the critical behavior of phase transition moves to smooth cross 
over. For degenerate quarks (m^ = rrid = nis = 0), the phase transition is again first order but 
with Tc ~ 155 MeV. For light up and down quarks and massive strange quark, the transition is 
cross over and Tc ~ 170 MeV again. Therefore, we would expect an upper value for the quark 
mass to secure first order phase transition. Above this values the order of phase transition seems 
to be cross over or weak second order. 



2 



It has been observed that the chiral symmetry is restored at the same critical temperature 
Tc ~ 154 — 174 MeV (depending on the quark flavors) as that of the deconfinement phase 
transition. The restoration of chiral symmetry means that the effective mass of quarks forming 
the confined hadronic states becomes zero. Another important consequence of chiral symmetry 
breaking restoration is the disappearance of mass degeneracy of hadronic states having same 
spin but different parity quantum numbers. 

We find that the bulk thermodynamic quantities at very hight temperatures of 4 — 5Tc 
remains below the Boltzmann limit jl51 114j : esb ~ gir'^T^/SO. This behavior would indicate 
that thermodynamic quantities may remain constant at much higher temperatures. This means 
that the deconfined matter might remain strongly correlated. That QGP at 4 — 5Tc turns to 
be strongly correlated, would have consequences on its hydrodynamic evolution, production of 
confined hadrons and might delay the freeze out processes. 

The chemical freeze out is characterized by s/T^ = 7 for three quark flavors [1]. s is the 
entropy density. This condition characterizes a stage at which annihilation and production 
processes are in chemical equilibrium. As shown earlier, the consequence that Tfo is below Tc, 
would lead to steep T-dependence of the inelastic scattering rate. 



4 Cosmological Consequences 

4.1 Nature and Order of QCD Phase Transition(s) 

At temperatures higher than the QCD critical temperature, Tc, the matter are mainly formed 
in quarks, gluons, leptons and photons. In QCD epoch, the energy scale {Aqcd ~ 200 MeV) 
is much larger than the physical masses of these components, we therefore can approximate 
their masses to be almost vanishing. As given above, the universe in this epoch has a radius of 
~ 10 km, which leads to mass content of ~ 4ttR%€{Tc)/3 ~ 1.25Mq. e{Tc) is the energy density 
at the critical temperature Tc- Its value is taken from lattice QCD simulations [13]. Under these 
circumstances the matter can be treated as a radiation. 

The relaxation time scale of particle interaction at QCD energy scale (F^ ~ a^T and Tg ~ 
—agThig, where T = g and Ua being gauge and strong coupling constant, respectively) is 
much shorter than the Hubble radius H, thus the different phases of the matter; QGP, hadron 
gas, leptons and photons, are likely in thermal and chemical equilibrium. Therefore this matter 
is much similar to radiation fluid and the effective d.o.f are the baryon quantum numbers. Under 
these conditions, the baryon number density can be calculated from 

nB{Tc) = V ■ s{Tc) ^ (1) 

S BBN 

where is the photon density at Big Bang Nucleosynthesis (BBN). r/ is the celebrated ratio of 
baryon density asymmetry (n^ — ng) to photon density. According to recent WAMP data [16j, 
rj reads 

4 X 10"^° < r/ < 7 X 10"^° (2) 

Using recent lattice QCD results [6j for Tc and s, we can calculate the physical units of baryon 
number density nB{Tc) for three dynamical quark flavors at physical masses in QCD epoch 

nB(Tc) = 8.1 X 10"" s(Tc) = 5.5 x 10"^ MeF^ (3) 

To get this value, we assume that no annihilation process is effective. The number of baryons 
filling out the Hubble volume up to the causal horizon at the QCD phase transition reads 

A^b(Tc)|^ ^ 3x10^8 (4) 
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Obviously, this number strongly depends on nature and order of QCD phase transition as s and 
Tc do (review Sec. [3j). When the universe loses ~ 30 MeV temperature, i.e. at the chemical 
equilibrium freeze out [1], Nb{Tc) consequently decreases to one half or one third of this value. 

In first order phase transition, the two phase co-exist and bubbles "dirt objects" in form 
of mixed phase of QGP and hadron gas are quite likely expected. To keep the temperature 
constant, bubbles slowly dominate the system. Since the two phases have different entropy 
densities, the bubbles release latent heat and accordingly the energy density e decreases |10) 
until the transition is completed. Released latent heat can reheat the matter (universe) to tiny 
fraction. Candidates for such "cosmic" bubbles would be primordial monopoles, cosmic strings, 
black holes, etc. As there is so far no verification for any of these objects, we have to look for 
other consequences, like primordial temperature fluctuations [17J and inhomogeneities forming 
dark matter clumps, etc. 

At QCD phase transition, neutrinos move freely with mean free path 
^/ree ~ 10^^ Rh ~ 1 cm and the fluctuations on their diffusion scale might be washed 
out during QCD phase transition, Xdiff ~ W~'^Rh ~ 1 m. These scales determine whether 
the temperature fluctuations (inhomogeneities) have to be neglected and the bubble hadroniza- 
tion proceeds in a homogeneous rate or the bubble hadronization is indeed inhomogeneous. 
According to [17J, the second scenario can be fulfilled, if the rms temperature fluctuations 
> 10~^. In that case, the scale of inhomogeneity in baryon distribution [T8] can be related 
to the scale of inhomogeneity in radiation fluid after the QCD phase transition. This might 
modify the neutron to proton ratio and thus explain the initial conditions for BBN. BBN is 
an important tool to verify our ideas about the early universe. Inhomogeneities would lead 
to formation of cold dark matter clumps in connection with the Hubble scale as given above. 
The scale separation between bubble hadronization and hadron diffusion is about two orders of 
magnitude. This scale difference has been calculated for first order phase transition. It depends 
on released latent heat and free energy difference between the two phases (QGP and hadron 
gas). 

Once again, in the first order phase transition, both phases exist at one critical pressure 
and therefore the speed of sound = dp/de vanishes throughout the whole transition. During 
the time of transition, the density perturbation within visible horizon (A < Rh) falls freely, i.e. 
radiation fluid velocity remains constant. At T > Tc, the density perturbation modes causes 
acoustic oscillations. They entirely vanish, when Cg vanishes. At very high energies, Cg = l/\/3 
in radiation-dominant matter (universe) |19j . Fig. [T] shows results for p/e and s/n in a hadron 
resonance gas. p/e can be related to Cs at very high temperatures, i.e., in an ideal gas. We 
notice that p/e reaches 0.14 at the critical point. 

The production of sterile neutrinos occurs through collisions between active neutrino gas 
and plasma of weakly interacting particles. The production in cold, warm and hot dark matter 
scenarios [20] is maximum around ~ 130 MeV. Intuitively, one can think about a close rela- 
tion between sterile neutrinos and QCD phase transition. According to [2T], first order phase 
transition significantly enhance sterile neutrino relic densities relative to crossover one. 

QCD predicts another phase transition at almost same Tc of confinement-deconfinement 
phase transition; the restoration of chiral breaking symmetry, i.e., anisotropic rotation of right 
and left handed quarks. The breaking chiral symmetry results in the pions as the lightest Gold- 
stone bosons. Therefore, strong fluctuations on the pion fields known as DCC, disordered chiral 
condensates [22], are quite likely through the transition. If chiral phase transition takes place 
through an out-of-equilibrium process, it would provide a crucial tool to explain the primordial 
magnetic fields [23]. The origin of the primordial magnetic fields shall be studied by the Square 
Kilometer Array (SKA) j24]. The particle production can be studied in hadron resonance gas 
(below Tc). Left panel of Fig, [2] shows ratios of various particles as a function of ^/s |25j . 

The dynamical fiuctuations associated with strong first order of phase transition are likely 



4 




very large. The continues second order or cross over phase transition might wish out large part 
of dynamical fluctuations in the final state and do not provide out-of-equilibrium. On the other 
hand, the dynamical fluctuations are conjectured to slow down near the second order phase 
transition. This has been confirmed in classical systems, solid state physics. In quantum field 
theory, the long-wavelength (spinodal) modes will be quenched through the second order phase 
transition |26j . 

In right panel of Fig, O we notice that the dynamical fiuctuations smoothly increase with 
the center-of-mass energy ^/s. This agrees with lattice simulations for Ug fluctuations [27] and 
might support the conclusions that any anomalous phenomenon associated with non-equilibrium 
phase change likely would be washed out, if the phase change is smooth and does not cause out- 
of-equilibrium. 

In second order phase transition, the correlation length ^ likely diverges at the critical point. 
This divergence can be characterized by a critical exponent v. v has a universal value, ~ 2/3. 
It does not depend on the microscopic details of the system in transition. It only depends on 
the universality class. 

m = c(r,)(^^) ' (5) 

The relaxation time r also diverges at the critical point, r gives the time needed for a small 
perturbation to be in-equilibrium. 



r 




Assuming that the system cools down with a constant rate so that 

m = (7) 
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Fig. 2: Left panel: ratios of different particles as functions of ^/s [25]. The largest ^/s is related 
to early universe, where the net baryon number is almost zero. Right panel: the dynamical 
fluctuations of K/tt particle ratios ^/s [28]. The points give the experimental data from various 
heavy ion collisions. The curves show the results from a hadron resonance gas [10] . 



where tq is the quench time. From Eq. ^ and ([7|), the correlation length and relaxation 
time can be related to t/tq 



m = Co(:;^) ' (8) 

The fluctuations that survive the phase transition are depending on the freeze out time. The 
fluctuations are characterizing these fluctuations. The freeze out time can be obtained by solving 
the equation r(i) = t 

\ i/(i+m) 

ror^) (10) 

As shown above, the hydrodynamic evolution of viscous QGP fluid is slower than the ideal 
QGP, i.e. r(t) > t or any out-of-equilibrium process is slow and therefore can survive the 
phase transition. That we do not observe dynamical fluctuations would indicate that the phase 
transition is continuous (cross over). Fig. [2j 



4.2 Chemical freeze out 

At chemical equilibrium, s/T^ = 7 for three quark flavors and the freeze out temperature 
Tfo ~ 174 MeV. s/T^ = 7 is a universal condition describing all experimental data from heavy 
ion collisions at a wide range on incident energies [l]. Then the entropy in QCD epoch reads 

S = V -T^ ^ 2.05 • 10^^ (11) 

This number reflects that about ten orders of magnitude an increase in produced particles as a 
reason of QCD phase transition are expected. The multiplicities of produced particles in heavy 
ion collisions should manifest this increase. Should we take into consideration that the Universe 
meanwhile expands this number becomes larger. 
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4.3 Non-ideal QGP (strongly correlated QGP) 



According to [29], the hydrodynamic evolution of QGP fluid with dissipation due to shear vis- 
cosity has a slower rate than the ideal QGP fluid. The transverse expansion in the viscous 
fluid is stronger and faster than the one in an ideal fluid. The pion production is considerably 
enhanced in viscous fluid, so that larger viscosity straightforwardly leads to more pions. The 
particle production increases when the freeze out surface is extended and the distribution func- 
tion become non-equilibrium. The freeze out surface [30] can be parameterized and fitted to the 
experimental data. With extension we mean change. The conditions deriving the freeze out as 
a function of the incident energies are discussed in [1]. 

The dissipation has another crucial consequence. It might modify the elliptic flow. The 
elliptic flow reduces with the viscosity. The reduction tends to reach a saturated value at high 
energies (large transverse momentum). The ideal fluid shows almost opposite behavior. The 
elliptic flow increases with the transverse momentum. 



5 Summary and Conclusion 

Lattice QCD calculations with dynamical quarks and physical masses show that the phase 
transition at very small net baryon density, which is related to early universe is likely cross 
over or very weak second order. Phenomenological indications for such a continuous transition 
is quite likely an ambiguous task. As discussed earlier, the order and critical temperature of 
deconflnement-confinement phase transition(s) calculated in lattice QCD strongly depend on 
the number and masses of quark flavors. The order varies from strong flrst order to very weak 
second order or cross over. Consequently, the critical temperature takes values between 150 MeV 
and 200 MeV. The cosmological consequences strongly depend on the QCD phase transition(s). 

We used to assume that QGP can be treated as a free ideal gas. The recent RHIC results 
suggested that QGP is a fluid rather than an ideal gas. Also recent lattice QCD simulations 
support such a conclusion. Taking into consideration that quark-gluon plasma is strongly cor- 
related (sQGP) and has shear viscosity would lead to the consequence that its hydrodynamic 
evolution is slower relative to QGP with zero viscosity. Also the transverse expansion in such 
a fluid is much stronger than the one in an ideal QGP fluid. Phenomenologically, the particle 
production is considerably enhanced in sQGP. All these results would to some extent modify 
our ideas about the time evolution in early universe, for example. 
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